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I. INTRODUCTION 

Quantum electrodynamics in (l+l)-dimensions, the 
Schwinger model [IJ , continues to play an important role 
as testing ground for field theory ideas. In this paper we 
use the Schwinger model to study the quenched approx- 
imation for the chiral condensate. The quenched chiral 
condensate has long been believed, from various indirect 
arguments, to be an ill-defined quantity in gauge theo- 
ries in any number of dimensions. Calculations in the 
Schwinger model are much easier than in higher dimen- 
sional gauge theories, and this is our motivation for the 
present study. One must keep in mind, however, the fact 
that spontaneous chiral symmetry breaking is prohibited 
in two dimensions by the Coleman-Mermin- Wagner theo- 
rem ■ For this reason we also consider the unquenched 
theory, where chiral symmetry is broken explicitly by the 
anomaly so that the chiral condensate should be well de- 
fined. 

The quenched approximation, in which the fermion de- 
terminant is discarded when generating the gauge field 
configurations, was first discussed analytically in this the- 
ory by van den Doel 0. Some of the subtleties were 
subsequently discussed in Refs. and H As we shall 
review in the next section, the indications of disease 
in the quenched Schwinger model very much resemble 
those in higher-dimensional quenched theories, analyzed 
by means of quenched chiral perturbation theory 0, • 
This may not be surprising, since the bosonized form of 



the Schwinger model 8] is a two-dimensional analogue of 
a chiral Lagrangian. The trouble with the quenched chi- 
ral condensate then stems, in both contexts, from the fa- 
mous double pole in the singlet correlation functi on pi HI. 
When analyzed in the finite- volume e- regime [Tollll|. the 
quenched chiral condensate is seen to be plagued with a 
"quenched finite-volume logarithm" at one-loop order in 
four dimensions ^2- Taken at face value, that is, if one 
were to push the expansion beyond its region of validity, 
this could indicate a divergent condensate. The analy- 
sis of the quenched chiral condensate in the e-regime has 
guided us in our present finite- volume calculations. 

There have been many Monte Carlo calculations of the 
chiral condensate in the Schwinger model, both quenched 
and unquenched (see [13, E3, [lH EH E3, EH Ea, E3l for 
recent work) . The first to study the quenched condensate 
via the distribution of the lowest Dirac operator eigenval- 
ues were Farchioni et al. ^3- They found quite an odd 
result — agreement with one universality class of sponta- 
neous chiral symmetry breaking at small volumes, and 
with another at larger volumes. We shall return to this 
issue in detail below. Kiskis and Narayanan ^3 recon- 
sidered the problem recently and found evidence for a 
divergent quenched chiral condensate, with Dirac eigen- 
values that did not appear to fit any of the three possi- 
ble chiral symmetry breaking classes. We shall see how 
these last two papers can be reconciled. In the process 
we shall consider volumes far exceeding what has been 
studied earlier. 

Our paper is organized as follows. In the next section 
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we briefly review the issues surrounding the quenched 
chiral condensate, particularly in this two-dimensional 
setting, and show how analytical arguments favor an 
ill-defined, divergent quantity. In Sec. Ill we turn to 
our Monte Carlo simulations of the quenched theory. 
We compare numerical results for the distribution of 
Dirac operator eigenvalues with distributions based on 
the three possible classes of spontaneous chiral symme- 
try breaking. We show how the distributions drift as the 
volume is changed, and how this explains the results of 
Ref. Moreover, we find that the distributions appear 
not to converge to any fixed limit, indicating that the 
spectral density p(X) does not attain a finite value at 
A = 0. By way of contrast, we present in Sec. IV numeri- 
cal results for the unquenched theory (Nf = 1), obtained 
by reweighting the path integral with the Dirac determi- 
nant. Our results there are consistent with the detailed 
predictions of the e- regime based on the explicit breaking 
of chiral symmetry due to the anomaly Section V 

contains a brief summary of our results. 



II. THE QUENCHED SCHWINGER MODEL 

Analytically, one can treat quenching by means of 
the replica method that considers Nf identical fermion 
copies, and then sends Nf — > at the end. This can be 
done trivially in perturbation theory at the fundamental 
level; in more than two dimensions it can be done in chi- 
ral perturbation theory |7j . Because the fermion determi- 
nant is exactly calculable in two dimensions, it can also 
be done beyond perturbation theory in the Schwingcr 
model. This was first realized by van den Doel H, and 
we shall here briefly review his calculation (see also 0,0)- 
At the same time we shall make contact with the very 
similar calculation in four-dimensions, in quenched chiral 
perturbation theory based on the replica formulation. 

In the continuum theory, we consider the Lagrangian 



C = —F flv F^-Y,M 



i$ + m — g4)ipi 



47T 



(1) 

with gauge group U(l), coupling g, and Nf species of 
fcrmions. A two-dimensional #-term has been included as 
well. Since the fermion determinant is exactly calculable 
in two dimensions, the theory is in large measure soluble. 
A convenient representation of the model is its bosonized 
form Q where the Lagrangian density, upon an exact 
integration over the gauge potential, takes the form 

Nf 1 o 2 ( Nl " 2 



N f 

+cm 2 y^A/"cos(2\/7r0j). 
j=i 



(2) 



Here Af denotes normal ordering and c = e 7 /2-7r where 7 
is Euler's constant. 



Although there is no spontaneous chiral symmetry 
breaking here, the bosonized action J2J bears a strong re- 
semblance to chiral Lagrangians in four dimensions, with 
/ = 1/ ypii playing the role of a (dimcnsionless) pion de- 
cay constant |2l| . This becomes particularly clear when 
we consider the theory in an expansion around static 
(zero-momentum) modes. If we define the partition func- 
tion in a sector of fixed topological charge v by means 
of 11] 
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'2(9), 



(3) 



the integral over can be performed exactly. Upon defin- 
ing the Nf x Nf matrix U = diag{exp(i2y / 7r0 J )}, the 
terms that survive in the static limit yield a Boltzmann 
factor 



(detC/fexp ~-Vm 2 AfTr (U + U f ) 



in analogy with the result in 4 dimensions ^lj (here V 
denotes the finite two-dimensional volume). Note, how- 
ever, that the mass term is proportional to m 2 rather 
than to m as in four dimensions. This is crucial for un- 
derstanding the difference with respect to spontaneous 
breaking of chiral symmetry. 

The analogy to a four-dimensional chiral Lagrangian 
holds to any order in chiral perturbation theory. In par- 
ticular it is useful for understanding the quenched limit 
of the theory. The mass term and its normal-ordering 
prescription complicate matters slightly, and it suffices 
to consider the massless limit. In that limit we read off 
the diagonal scalar propagator from Eq. J3J), 



S(P 2 ) = 3 



p 2 7r p 2 (p 2 + N fg 2 jii) ' 



(4) 



which reduces to the classic Schwinger result of an or- 
dinary massive propagator when Nf = 1 0]; in the 
quenched limit (Nf = 0) a double pole develops, 

as first noted by van den Doel 0. This phenomenon is 
completely analogous to what happens in quenched chiral 
perturbation theory in four dimensions when formulated 
in terms of replicas; see the appropriate replica Feynman 
rule in Eq. (7) of Ref. 0- The only difference is the re- 
placement of the Schwinger mass parameter fi 2 = g 2 /it 
by what is there commonly normalized as /N c , where 
N c is the number of colors. (The double pole had of 
course been observed much earlier in quenched chiral 
perturbation theory by means of a supersymmetric ex- 
tension 0,0-) 

The above description of the diagonal scalar propaga- 
tor glosses over the fact that in two dimensions the propa- 
gation of massless degrees of freedom require an infrared 
regularization. This is true even when the Schwinger 
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mass fx = g / \pn is taken into account because of the re- 
maining l/p 2 -poles in the propagator Q when Nf — > 0. 
After regularizing this infrared divergence by an addi- 
tional mass parameter mm, the calculation of the chiral 
condensate requires integration over a closed loop of the 
propagator 



where 



G(p 2 ) 



p 2 + (N f - l)fi 2 



{p 2 +mj B )(p 2 + N f Lt 2 + mj R ) p 2 + M 2 



where M is the arbitrary mass defining the normal- 
ordering prescription ||. From here the quenched chiral 
condensate has been computed 0, 0, HJ to give 



(tpip) = - lim cM {Af cos [2 v^fc]) 



hm cmjB 1 H % — 

Nf->0 v iJ < 



1/2 JV/ 



(7) 



For fixed infrared cutoff m/R, the limit iV/ — > yields 3 

(^V> = - lim cm IR e>* 2/2miR , (8) 

which is infrared divergent. This ordering of limits seems 
closest to the actual "physical" (i.e., computational) def- 
inition of the quenched theory, and we shall view it as the 
simplest manifestation of the difficulty with defining the 
quenched Schwinger model. Since the result is divergent, 
and since the computation has been done based on an 
expansion around the massless theory, one can question 
to what extent the resulting divergence is a truly reliable 
prediction. More detailed computations at fixed V can 
be found in Refs.0 andl22l 

III. MONTE CARLO ANALYSIS 

We consider a lattice with volume (La) 2 and employ a 
non-compact formulation for the gauge field. In a sector 
with topological charge Q, the gauge field may be 
decomposed [23| as 



Ai(x) = d*(j>{x) + 



2* 
L 



di<x(x), 



(9) 



A 2 (x) = -d* x 4>(x) + ^h 2 - d 2 a(x) - (10) 

Here <fi(x) is a real periodic function with no zero mode, 
and /i M are two real constants in the interval (—1/2, 1/2] 
that parametrize the two Polyakov loops on the 2-d lat- 
tice. (9 M and d* are forward and backward finite differ- 
ences; a(x) represents the gauge degree of freedom. The 
gauge action is then 



(11) 



A0(x) = Y, W + M) + <K* - M) - 20(x)] . (13) 



We denote the gauge coupling by (3 — l/(2g 2 a 2 ). The 
continuum limit may be taken at fixed finite g and at 
fixed volume by taking (3 and L to infinity. Alternatively, 
we can keep g and L fixed and vary the physical volume 
by changing j3. 

Our choice of the non-compact action is motivated 
mainly by the ease with which we can generate inde- 
pendent gauge configurations with the desired topologi- 
cal charge. Successive configurations are generated by a 
heat bath so that there is no autocorrelation. We restrict 
ourselves to the sector with zero topological charge. 

We employ lattices of linear size L ranging from 8 to 60, 
and fix (3 = 2. We also have one data set with (3 — 1 
and L = 48; if we renormalize at fixed g as discussed 
above then this is equivalent to L — 48\/2 ~ 68 at (3 = 2, 
allowing us a larger physical volume at modest additional 
cost. 

We use the massless overlap-Dirac operator |24f for the 
fermions, 



(14) 



Here Hw = 75-Dw( — 1) is the hermitian Wilson-Dirac 
operator with mass parameter set to — 1. With this nor- 
malization the eigenvalues of the hermitian overlap-Dirac 
operator H = 75D lie in the interval [—1,1], but there 
is a wave-function renormalization = 2 with respect 
to the conventional normalization of the Dirac operator 
[25| that one needs to keep in mind. 

For L < 32 we diagonalize Hw exactly by a House- 
holder transformation followed by QL iteration [26| . 
From this we construct H, whose eigenvalues we obtain 
similarly. For the larger lattices we compute the lowest 
eigenvalues of H 2 with the Ritz variational algorithm of 
Kalkreuter and Simma |27|. Here the action of e(Hw) 
on a vector is obtained using a modified version of the 
two- pass algorithm In the first pass, the Lanczos al- 
gorithm is used to obtain a tridiagonal matrix, Tyy, that 
is a good approximation to Hw- Then e(TV) is obtained 
by an exact diagonalization of T\y, and the second pass 
is used to compute the action of e(Tw) on a vector. 

Previous papers have compared the eigenvalue dis- 
tributions to the predictions of random matrix theory 
(RMT), and we shall proceed to do the same. In higher 
dimensions, where spontaneous chiral symmetry break- 
ing is allowed, it is by now well known that the lowest 
eigenvalues of the Dirac operator in the e-regime are dis- 
tributed according to universal finite-size scaling distri- 
butions that can be derived either from RMT I291I3CI o r 
directly from the chiral Lagrangian framework |3ll l32| . 
There is ample evidence that the universality class of 
these distributions is dictated by the way the fermions 
transform under the gauge group. This has been demon- 
strated for three different gauge group representations 
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with overlap fermions |33| . and for a variety of exotic 
representations with staggered fermions [34| . 

According to this classification, if chiral symmetry 
could be spontaneously broken in the Schwinger model, 
where the fermions transform as a complex representa- 
tion, the relevant universality class in RMT terms would 
be that of the chiral unitary ensemble (chUE). It is by 
no means clear, however, whether these predictions are 
of any relevance to the quenched Schwinger model, where 
spontaneous chiral-symmetry breaking is prohibited 0- 
Ref. included such a comparison, but the results were 
not easy to understand. For small lattices the eigenvalue 
distributions appear to fall in with the chiral symplcc- 
tic ensemble (chSE) , which should be of relevance to real 
fermion representations. For larger lattices there appears 
to be a switch to the perhaps more natural chUE. With 
our greater statistics and our larger lattices, we are now 



able to sec that such surprising results do not really hold. 

We begin with attempts to fit the distribution of the 
lowest eigenvalue to the form predicted by RMT for the 
quenched (N f = 0) chUE 

Pi«i) = ^ie- £?/4 , (15) 

where £i = AiEL 2 and the condensate E is the fit param- 
eter. We compare to histograms with 20 bins in £i. The 
comparisons are shown in Fig. ^ for nine volumes, and 
the results of the fits to the RMT prediction are listed in 
Table 13 

L = 24 is the case that comes closest to agreement with 
the chUE prediction, but even here our fit gives a x 2 /dof 
of 58.1/19 corresponding to a confidence level of 7.6 x 
1CP 6 . Our high-statistics data enable us to rule out the 
chUE scenario here. As can be seen from Table [U for all 
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TABLE I: Number of configurations studied for each volume 
V — L 2 , the condensate E, the \ 2 1 number of degrees of 
freedom, and confidence level from fits of the lowest eigenvalue 
distribution to the RMT form of the chUE with 20 histogram 
bins. The last line is for /3 = 1, all others for /3 = 2. 
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other lattice sizes the RMT fits to the chUE predictions 
are ruled out even more thoroughly. 

Farchioni et al. found agreement with the chUE 
prediction for large volumes, in particular for L = 16 and 
(3=1 (see their Fig. 6; their definition of (3 is twice ours). 
This is the same physical volume as L ~ 23 and (3 = 2. 
As we stated above, we rule out the chUE at this volume 
and at all other volumes studied. 1 For smaller volumes, 
Farchioni et al. favor the chSE, claiming in particular a 
good fit at L — 16 for (3 = 2. A fit of our high-statistics 
L = 16 data to the quenched chSE distribution gives 
^ 2 /dof = 436/19, ruling it out as well. 

As seen in Fig. ^ the peak of the distribution of the 
lowest eigenvalue moves downward faster than 1/L 2 , the 
RMT prediction, while a sizable tail of the distribution 
persists. The RMT scaling law is based on the assump- 
tion of a finite eigenvalue density p(0) at the origin; a 
scaling faster than 1/L 2 thus entails a divergence in p(0) 
as L — ► oo . By the Banks-Casher relation £ = 7rp(0), 
this in turn implies a divergent chiral condensate. In- 
deed, from Table [I] we see that the fitted £ grows quite 
rapidly with increasing L. 

We now elaborate on this last point. Let Xi(L,(3) be 
the i-ih lowest non-zero eigenvalue of H on an L 2 lattice 
at a coupling (3. For a finite chiral condensate to form 
we expect that 

approaches an i-dependent constant as L — > oo, 

lim ML, (3) = M(3) . (17) 



Ref. 1 1 El describes a simulation of a compact gauge action, how- 
ever, so the comparison between our calculations cannot be made 
precise. 



In lattice units we expect f t ((3) oc £ |H|. 

Furthermore, 

the products fi((3)V(3 should approach finite continuum 
limits for (3 — ► oo. 

The scaled variables fi(L, (3)\/j3~ and fc(L,(5)\[(3 are 
plotted as functions of the physical size L/y/]3 in Fig. El 
(We combine data for (3=1. with data for (3 = 2.) We 
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FIG. 2: Evidence for a diverging condensate in the quenched 
Schwinger model. 



see that these quantities do not approach L-independent 
constants as one would expect on the basis of the exis- 
tence of a finite chiral condensate £. 

The distribution ^2(62) of the second scaled eigenvalue 
(with £ taken from Table P) for some large volumes are 
shown in Fig. [3] and compared with the predictions for 
the (Nf = 0) chUE [U , 

P 2 (6) = \&~ eji J q 2 duu [•*!(«) - h(u)h(u)] . (18) 

Much like the distributions of £1 above, these clearly do 
not fall in the universality class of the chUE. We can 
eliminate the scale £(£,/?) in these comparisons, by plot- 
ting the distribution of r = \i((3, L)/ \2((3,L). We com- 
pare this to the prediction from the chUE [3ol |3(j 

r , 1 r r j 3 ( u2 \ 

P\ r ) = 77; 777 / duu exp -— — 

4(l-r 2 ) 2 7 4(l-r 2 )y 

x [I 2 {u)-h (u)I 3 (u)] , (19) 

in Fig. 2| Again, the data do not fall in the universality 
class of the chUE. 

Do the Pi(£,i) reach some limiting distributions as L — > 
00? It is difficult to answer this question based on the 
data shown in Figs.nand|21 It is conceivable that 
approaches a function peaked at zero while P2 (£2 ) reaches 
a limiting form peaked away from zero, but we have no 
real evidence for this. Whatever the answer, however, 
it is possible that the distribution p(r) of the ratio does 
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FIG. 3: Distribution of the scaled second eigenvalue of the quenched Schwinger model. 
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FIG. 4: The distribution of the ratio of the two lowest eigenvalues in the quenched Schwinger model shows that they do not 
obey the universal distribution given by the chUE. 



have a limiting distribution. Even though Ai and A2 
go to zero much faster than 1/L 2 , we can discern some 
level repulsion that favors a ratio r between and 1. 
Figure shows (r) as a function of the physical size and 
this average seems to approach a finite limit as L — > 00. 



IV. THE UNQUENCHED THEORY 

When we compute all the eigenvalues of the Dirac oper- 
ator on each pure gauge configuration, we can calculate 
observables in the theory with Nf ^ by reweighting 
with the fermion determinant. We restrict ourselves to 
Nf = 1 since statistical fluctuations are worse when the 
target ensemble is farther from the original quenched en- 
semble. In the continuum Nf = 1 is of course Schwinger's 
original model Q], which is exactly soluble. In particular, 
the infinite-volume chiral condensate can be computed 
analytically, @ 
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= (0.1599...) 5 , (20) 

in the conventional normalization, where /1 = gj\pH is 
the Schwinger mass and 7 is Euler's constant. Because 
the chiral symmetry is broken explicitly, the analysis of 
Leutwyler and Smilga [ljj and the whole RMT analysis 
for the Nf = 1 theory apply directly here. We therefore 
know the complete microscopic spectrum of the Dirac 
operator in the e-regime and it belongs to the universality 
class of the chUE. 

Once again, the simplest quantity with which to com- 
pare is the distribution of the smallest (non-zero) Dirac 
eigenvalue [35|. 

= \txh{!ix)e-Z'\ (21) 

here restricted to the massless case. We were able to 
reweight ensembles on lattices up to size L — 32. Again, 
we made fits to the RMT prediction, with E the fit pa- 
rameter, using histograms with 20 bins. The fits are 
detailed in Table [H] and shown in Fig. The agree- 
ment with the chUE is quite good for the two largest 
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20 30 
L/V0 

FIG. 5: Expectation value of r = X± (L, {3)/\ 2 (L, (3) in the 
quenched Schwinger model as a function of the physical size. 
The data are not consistent with chUE but do seem to ap- 
proach a finite limit as L — » oo. 




FIG. 6: Comparison of the distribution of the rescaled lowest 
eigenvalue for Nf — 1 with the prediction of RMT for the 
N f = 1 chUE. 



TABLE II: Number of configurations studied for each volume 
V — L 2 with ft = 2, reweighted to Nf = I, the condensate 
E, the x 2 ) number of degrees of freedom and confidence level 
from fits of the lowest eigenvalue distribution to the RMT 
form with 20 histogram bins. 
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lattices. The second eigenvalue (Fig. EJ) tells a similar 
story. We show the averaged ratio (r) in Fig. |S| which 
may be compared to Fig.[S]for the quenched theory. The 
result for L = 32 agrees with the prediction of the chUE, 
(r) = 0.5044 HI. 

Converting the result from Table [n] for the largest 
lattice size to the conventional normalization, we find 
Ti/g = 0.1648(11), within 3% of the continuum value 
given in Eq. Ij20|l : finite lattice spacing corrections, ex- 
pected to be of O(l/0), could easily explain the differ- 
ence (see also 0|). 2 We note, however, that fits to the 
RMT distribution of the lowest eigenvalue do not work as 
well for the smaller lattices, L < 24, giving unacceptably 
small confidence levels. It would be interesting to see 




FIG. 7: Rescaled second eigenvalue for Nf = 1, compared 
with the prediction of RMT for the Nf = 1 chUE [see Eq. (15) 
inRef.il]. 



2 The fit to for L = 32 gives a value of E that is 7% larger 

than the value shown in Table ITT1 This is consistent with wide 
experience that higher eigenvalues suffer larger finite-volume ef- 
fects. 



whether the agreement with RMT persists on larger lat- 
tices and whether the continuum value for S is correctly 
reproduced as (3 — > 00. Unfortunately, reweighting be- 
comes prohibitive for larger volumes. A direct numerical 
simulation of the Nf = 1 theory is probably needed. 
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FIG. 8: Expectation value of r = Xi(L, (3) j\i (L, /3) for N f 
1 as a function of the physical size. 



V. SUMMARY 

As we have shown, the quenched Schwinger model does 
not fall into any of the three universality classes of chiral 
symmetry breaking. In view of the Coleman-Mermin- 
Wagner theorem, which forbids spontaneous breaking of 
a continuous symmetry in two dimensions, this in itself 
is not very surprising. The only possible loophole out 
of this argument would be to note that the quenched 
theory is non-unitary and thus it might not satisfy some 
assumptions of the theorem. 

The application of the Coleman-Mermin- Wagner the- 
orem to the quenched theory may be considered in two 
ways. In the so-called supersymmetric formulation of 
quenching, spontaneous chiral symmetry breaking is as- 
sociated with quenched Goldstone bosons and Goldstone 
fcrmions. In the replica formulation, the spontaneous 
breaking is associated with Goldstone bosons alone. In 
both ways of considering quenching, the theorem seems 
to exclude rigorously the possibility of a non-zero con- 
densate. In that light it is perhaps surprising that the 
numerical evidence now points towards an ill-defined, di- 
vergent chiral condensate, rather than a vanishing con- 
densate. 



The Schwinger model with Nf = 1 is on an entirely 
different footing due to the explicit breaking of chiral 
symmetry by the anomaly. Here we have unambigu- 
ous analytical predictions for the behavior of Dirac op- 
erator spectra near the origin, and our Nf = results 
re-weighted with the Dirac determinant to simulate the 
Nf = 1 theory are consistent with these analytical pre- 
dictions. For the massive Schwinger model with Nf = 1, 
a disagreement with RMT must appear as the mass is 
taken to infinity. Our statistics have not been good 
enough to attempt the even more ambitious reweight- 
ing to simulate the Nf — 2 theory. Also here unusual 
results should appear, presumably with the Dirac oper- 
ator eigenvalues being strongly repelled by the origin so 
as to produce a vanishing p(0). For Nt = 2 analytical 
calculations 01 suggest a behavior 



p(A) ~ AV3 



(22) 



and an interesting question is whether such behavior at 
the rescaled level has a universal distribution from "criti- 
cal" Random Matrix Theories (the precise behavior 11' I'D 
is actually realized in a very simple chiral matrix model 
or analogous eigenvalue models |38||. This could be 
an interesting topic for future investigations. 
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